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La Comète Orkisz (1925 C). 


Sa découverte. Son orbite. 
Formules de nouveau genre appliquées à la détermination d'une orbite. 
Par Thadée Banachiewicz. 


La Comète fut découverte à la Station astronomique provisoire, située au mont Łysina (h = 912 m, 
À = 20°4, ® = 49"46) de la chaîne du Beskid, aux environs de Cracovie. Le programme d'obser- 
vation de cette station comprend l’investigation des conditions atmosphériques au point de vue des exi- 
gences de l'astronomie, les observations d’occultations d'étoiles par la Lune, les étoiles variables, la recherche 
des comètes. 

La découverte fut faite à l’aube du 3 Avril 1925 par M. Lucien Orkisz, assistant de lOb- 
servatoire de Cracovie. Il ne fut guère possible de déceler ce matin-même un mouvement propre, mais M. 
Orkisz constata au Catalogue de Dreyer l’absence d’une nébuleuse pour ce point du ciel. 

Vu la communication assez difficile avec la Station, M. Orkisz tâcha le 3 Avril de se relier 
à l'Observatoire de Cracovie à l'aide d'un héliotrope, mais ce fut en vain. Le 4 avril il constata que l'objet 
suspect avait changé de position ce qui présentait une confirmation définitive de son caractère cométaire. 
Examinée au refracteur Merz (objectif de 115 mm, gr. 46) la Comète se présentait comme une nébuleuse 
presque circulaire d'un diamètre de plusieurs minutes. Vers 3° 13° T. U. la Comète était déjà à peine visible. 

L'Observatoire de Cracovie obtint la nouvelle de la découverte le soir du 4 Avril et la transmit 
aussitôt télégraphiquement à la Centrale de Copenhague. - 

La position de la Comète communiquée à la Centrale: x — 22" 27", ò — + 16° 14° présente 
cependant non la position observée, mais celle prévue par le découvreur pour le 5 Avril d’après les lectu- 
res des cercles de l'instrument, faites pendant les observations du 3 et 4 Avril; la brièveté de la dépêche 
de £ysina était la cause de l'omission de ce détail. 

Le 5 Avril M. K. Kordylewski, s'étant rendu à la Station de Lysina, détermina d’après les 

croquis de M. Orkisz les positions suivantes de la Comète. 


PO RO 20 0. ango —— 22 22" 14e Osip — T 1423 .2- (découverte) 
FAUNE . » :-$ : 22 23 40 y +15 21.9 


M. Kordylewski estime pour les limites des erreurs de ces coordonnées +8° respectivement? 2”. 

Le 4/5 Avril le ciel était couvert à Cracovie, mais les observations de la Comète ont pu être 
faites en Pologne à Varsovie (M. F. Kępiński) et à Wilno (M. S. Szeligowski et M Kowal- 
czewski), le lendemain à Poznań (M. S. Andruszewski) et Léopol (M. L. Grabowski). 

Les observations ci-dessus de M. Orkisz, combinées aux autres, ont servi à M. J. Witkowski 
à calculer une orbite approximative, sur laquelle était basée notre première éphéméride, télégraphiée à Co- 
penhague le 11 Avril. 

Détermination de l'orbite et calcul de l'éphéméride. 


Nous allons exposer ici avec quelques détails notre calcul de l'orbite et les formules qui nous ont 
servi, parce qu’on ne trouve pas dans la littérature du problème un ensemble de formules arithmo- 
métriques, et parce que nous croyons utile de donner aux calculateurs un exemple numérique de l’appli- 
cation des formules astronomiques de nouveau genre. Nous rappelerons d’ailleurs aussi un 
procédé simple pour trouver la valeur approximative de la distance d'une comète. 

Dans la détermination de l'orbite nous nous basons sur les observations suivantes 


Temps d'observation mn PEUC: Diop gÉOC. 
No. 1 1925 Avril 5.1161 T. U. 22"26" 43.30 <+16°37/20”.0 Obs. F. Kepiñski (Varsovie) 
No. 2 ~ m 8.1138 , 29 42.66 +19 46 28 7 .„ J. Witkowski (Cracovie) 
No. 3 as PTIT UD T, 32 mpn AGa 23 4559 93 5 S 


La première observation nous a été aimablement communiquée par M. le prof. M. Kamieński, 
directeur de l'Observatoire de Varsovie; nous l'avons déjà corrigée de l'effet de la parallaxe (—0°.21 +4”.0) 
et de l'effet de l'aberration des étoiles fixes (+1°.02 4-11”.0). Les corrections correspondantes ont été 
introduites aussi dans les observations No. 2 et 3 de M. J. Witkowski. Pour tenir compte de l'aber- 
ration planétaire, on a appliqué dans le calcul de l'orbite les corrections —0°.00100, —0*.00099 et 
— 0100098 aux moments des observations, ces corrections ayant été obtenues d’après l’'éphéméride que 
nous a envoyée M. Kamieński. 

Voici donc définitivement les données que nous avons employées dans la détermination de l'orbite. 


t. d’obs. — t. d’aberr. Cosinus directeurs 1925.0 Coord. rectang. du Soleil 
a b c X Y Z 
No. 1 1925 Avril 5.1061 +0.87994 — 0.37932 <+0.28606 +0.96737 +0.23477 +0.10184 
No. À. ” 8.1039 +-0.86894 —0.36121 <+0.33832 <+0.95375 +0.28032 —+0.12160 
No. 3 f » _11.0991 +-0.85432 —0.34122 “<+0.39205 <+0.93763 +0.32509 +0.14102 
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On calcule tout d’abord 


di X & ay Xa à) 
Ps : Pi = M = — | bi Yy b| (ét) : | bs Ya dl (h-t). 
Eci Za C9 C3 Z2 C3 


On remarque que les deux déterminants ont deux dernières colonnes communes, on les déve- 


loppera donc suivant les termes de la première colonne. 
On trouve M = 108396 : 113808 = 0.95245. 


D'après les formules du type x = —X -|-ap on a, pour les coordonnées héliocentriques de la 
Comète à la première et à la troisième date, les expressions numériques suivantes 

x, = — 0.96737 +0.87994 p, y, = -0.23477 - 0.37932 ọ; z; = —0.10184 +0.28606 p, (1) 

x, = — 0.93763 +0.81370 p, y, = — 0.32509 — 0.32499 p, 23 = 0.14102 +0.37341 p, 


d'où l'on tire r, ? = 1.001293 — 1.582614 p, +1.000008 p, ? = R,? — 2 R, cos% p, +p,? (2) 


ra ° = 1.004720 — 1.419914 pọ, +0.907161 p, °? = R,” — 2 R, M cos h; pı +P:° 
(xs — x)? + (y, g y) + (z; mi O 0.010577 E 0.020599 Pi +0.014970 Pi 2 = E 2 hg COS @ Pr -+ h° Pi ï 


Valeur approchée de la distance. 


D'après la note dans les Comptes Rendus de l'Académie des Sciences de Paris, t. 161 (1915), 

p. 122, on a, pour la distance pọ de la comète à l’époque moyenne, une équation très approchée 
(p* — 2pacosp.+ a VR?—2Rcosbp + pọ? = 22? R?, (3) 
où lon a posé a=g{i+M):2h. 

Si l’on construit sur un plan le point R d'après les coordonnées rectangulaires a cos et asin 9, 
et le point S d’après les coordonnées À cos ® et Rsin®Ÿ, il suffit pour la détermination de pọ de trouver 
la distance de l'origine d’un point C, situé sur l’axe des abscisses, pour lequel 

CRÉES TR (4) 

En défaut d'un diagramme des courbes de Cassini (généralisées), ce problème géometrique se résout 
très rapidement par les mesures des longueurs CR et CS pour différents points C sur l'axe des abscisses. 
Le point pour lequel le produit CR? X CS aura la valeur donnée 2 4° R”, donnera l'abscisse cherchée p. 

Dans le cas numérique examiné on a d'après les coefficients de la dernière équation (2) 

g? = 0.010577 h*= 0.014970, par suite a = 0.8202, et comme 29ghcos® = 0.020599 (d'après la 
même équation), on obtient 
acos® = 0.672 a sin p = 0.470. 

D'autre part, les observations étant presque équidistantes, on peut identifier À avec K, et cos Ÿ 


avec cos Ù», de sorte que 
R cos Ÿ = a, X, + b, Y, + ©, Z, — 0.7686 (5) 


Comme R, = 1.002, on trouve cos Ÿ = 0.7671, d’où sin p = 0.6415, et FR cosp = 0.769 R sin d = 0.643. 
Dans le cas général on trouverait Rcos® comme moyenne arithmétique de À, cos, et R, cos h 
On a obtenu enfin 2a°kR* = 1.351. 
En prenant maintenant le point C sur l'axe des abscisses à distances de l'origine égales à 1.60, 
1.70 et 1.80, on obtient d’après les mesures sur le diagramme : 


pour ọ = 1.60 CR = 1.040 CS = 1.051 d'où le résidu de l'éq. (4) —0.214 


1.70 1.131 1.132 Mn IE Re aT 
1.80 1.224 l'0187 SRE Rd T O 
1 
En appliquant une fois la regula falsi on obtient ọ = 1.669. Comme, d'autre part, 6 = Que Pi, 
il suit p, — 1.710 (approximativement). 


Tout ce calcul ne demande que quelques minutes. 

D'après la construction il est évident qu’il n'existe aucun autre point C, satisfaisant à l'équation (4); 
la solution est donc unique. Ce résultat était d'ailleurs à prévoir parce que, d'après le théorème de notre 
note déjà citée, les équations de la méthode d'Olbers n'admettent qu’une seule solution, quand la distance 
angulaire de la comète au Soleil dépasse 19° et le mouvement apparent diurne sur la voûte céleste est 
supérieur à 8.5; c’est justement le cas qu'on rencontre le plus souvent en pratique. 


On arrive maintenant à la valeur exacte de p, en calculant S4 — S, d’après les formules usuelles 
pour p, = 1.70 et p, = 1.72. La condition s; — s, = 0 donne alors ọ, = 1.71202. On remarque que la 
valeur de p, fournie par la construction présente une bonne approximation. 

Avec la valeur définitive de p, on obtient d'après ( ) 


Le —= +0.539105 y = — 0884173 z; = +0.387900 

x, = +0.455441 Y, — — 0.881479 z, = +0.498265 
Dans le procédé suivant de la détermination de l'orbite, qui est bien préférable, en cas de calculs 
arithmométriques, aux méthodes des Traités, on calcule maintenant 9 quantités P,....... R,, qu'on pourrait 


nommer éléments vectoriels de l'orbite. On a 


R; = (Yiz — y; zı): S Ry = (zi xy — z; x) >S R, = (xi Ys ~ xX; Yı): S, (6) 
oS E + G a HTF (y p) 
Numériquement S = 0.153107 1 : S = 6.53138 
R, = —0.64417 k, = —0.60057 k, = — 0.47366 
L'angle 2f entre r; et r, s'obtient d’après les formules 
ERE I a a O e A U Us Tuez 


et comme r, r cos 2f = 1.21819 on trouve tg 2f= 0.125684, 2f= 709'48”.96, f= 38034'54”.48, 
sin f = 0.062473, cos f = 0.998047. 


; ps 1 
Pour calculer les anomalies v, et v, désignons v = 5 ( + va). 


On a évidemment v, = v — f, v, =v + f et on trouve 
r= -T | 
tg Le si iala ctg T ou mieux tgv = EURE PUS Se (7) 


+ Vr, =+- Vr, 
Si lon compte v de — 180° à 180°, v a le signe de rą — r, ce qui détermine son quadrant. 
D'après la dernière forınule (7) nous avons obtenu numériquement tg v — 27.74 : 432.8 — 0.06409, 
v — 3°40'41” (au lieu de 01”). Il en vient v, — 0°5/46/’.5, v, — 7°15'35”5 et, d’après la formule 
q = r cos? li v, q — 1.10583 (d'après v,), et q — 1.10582 (d’après v); en moyenne: 


q — 1.10582, 


2 ir rs — (ri + r3) cos f 


D’après v, et q et d’après v, et q on obtient par le procédé connu deux valeurs de Tt (l’époque 
du passage «u perihélie), qui doivent coincider. Numériquement on trouve tT — 1925 Avril 5.0258 et Avril 


5.0261, et, en moyenne t — 5.0260 


Envisageons les axes héliocentriques x et y dirigés suivant les rayons vecteurs v == 0° et 
v — 90°, et introduisons 
X = T, COS V, y =r sinv, X} =T} COSV} Ya = T} SİN V3 (8) 


ga: (xı Y3 — X yı). 
Ceła étant, la formule de nouveau genre 
RP E rer) A om EAA (9) 
|Q: Q, QÍ pu Lx y zf | \—0cy ox, J 


donne immédiatement P.....Q,; ce sont les cosinus des angles que forment les axes x y avec x y z. 
Remarque. La relation (3) pourrait être écrite sous forme de six équations ordinaires 


Ps = o (x, Ya — Xs yY) P, = o (Y, Ya Ua Hi). Pr = o (Z, Ys — 2s Yı) (9*) 


Q: = o (—x, x + X3 xX,) Q, = o (—Y, X3 + Y; xX) Q: = o (—2z, X3 + zx) 


Non seulement les équations (9*) sont plus compliquées à écrire, mais encore — et cela est le point essentiel— 
demandent-elles un effort mental notablement plus grand de la part du calculateur, que quand on fait usage de la for- 
mule (9). Les avantages des formules de nouveau genre sont d’ailleurs encore plus accentués, quand les transforma- 
tions à effectuer sont plus complexes. 


Numériquement on trouve 


|P- BR P, | | +0.539105 — 0.884173 +0.387900| ( +0.916385 —7.193511 0| [048850 — 0.79955 +0.34942] 
Q. Q, Qı} = 4 +0.455441 —0.881479 +0.498265, {—0.012135 +7.222641 0; = | —0.58857 —0.00630 +0.80843; 
RAR PA |—0.64417 —0.60057 —0.47366 | | 0 0 1 |—0.64417 — 0.60057 —0.47366| 


où l’on a ajouté la dernière ligne, calculée déjà ci-dessus, pour avoir tous les éléments vectoriels ensemble. 
Les éléments P,....k#,, avec q et t déterminent complètement l'orbite et suffisent pour résoudre 
tous les problèmes qui dépendent des éléments de l'orbite. En particulier, ils permettent de calculer les 


éléments ordinaires &, 6}, à d'après les formules 


R, = sin i sin Q, k = — sin Í cos Q, : 
aE. Mr E R, = cos i. (10) 
P, = sin į sin ®, Q, = sin į cos Ō, 


Si l’on désirait rapporter les éléments à l’écliptique, on le ferait par la formule 


(P: Q: LE p Q: R2) P, Q, £ | 1 0 0 | j + 0.48850 — 0.58857 — 0.64417 
P, Q, E} =P, Q, R,?. p (€) = p Q, R, $ 0 cose— sine ¿=< — 0.59447 + 0.31592 — 0.73945, (11) 
i Q. RT MPNO R, J P, Q, R?) (0 sine me) pa 0.63873 + 0.74417 — 019556] 
écliptique équateur équateur 
On a posé dans ce calcul sin € = 0.39793, cose — 0.91741 (1925.0); v. Supplemento Internationale 
Nr. 3, p. 51, du Rocznik Astronomiczny Obserwatorjum Krakowskiego. 
On obtiendrait maintenant aussitôt au moyen des formules (10) 


w = 40°38’ oe 315756. a T (1925.0) 
[l suffit d'écrire ces éléments traditionnels en nombres ronds, puisqu'on n'en a pas besoin ni 


dans le calcul de l’éphéméride, ni dans la correction de lľorbite, ni dans le calcul des perturbations. Les 
éléménts vectoriels les remplacent partout avantageusement dans ces opérations. 


On peut réduire les éléments vectoriels à une époque quelconque au moyen de la multiplication 
du cracovien des éléments vectoriels par le cracovien, dont on donne les termes dans le Supplemento 
Internationale du Rocznik Astron. Obserw. Krakowskiego. 


Les coordonnées rectangulaires héliocentriques s’obtiennent d’après la formule de nouveau genre, 
que nous écrivons ici sous sa forme générale 


EE fab} BE A 
fl =) 1I g aro (12) 


La formule (12) est évidemment bien plus commode pour le calcul arithmométrique que les for- 
mules classiques qui demandent la formation de A+v, B+v, C+v et la recherche des sinus de ce 
trois angles. En outre, les constantes de la formule (12) s'obtiennent au cours de la détermination de 
l'orbite, tandis que les constantes de Gauss demandent un calcul spécial, du moins dans les procédés les 
plus recommandés. Par la même raison le procédé exposé est plus commode que les formules arithmomé- 
triques en usage en Grande-Bretagne (voir The Journal of the British Astronomical Association, Vol. XXXII, 
Nos. 6 et 8), qui supposent la détermination préalable des éléments ordinaires, ©, $} et i; étant „une 
formule de nouveau genre“ la formule (12) est d’ailleurs aussi un peu plus compacte et plus facile 
à l'usage que les formules anciennes. 

En faisant abstraction de la forme nouvelle, lexpression (12) se rapproche encore le plus des 
formules de Gibbs (On the determination of elliptic orbits from three complete observations), mais ce 
géometre éminent n’a pas reconnu la formule simple (9) pour le calcul des constantes P,....0.. 

Pour une orbite parabolique on a 


(13) x=q (1 — m’) y — Ÿq m m +4 m°— M (l'équation de Barker). 
1 
Les tables pour obtenir m= tg 3 v d'après M sont seulement en préparation; en attendant qu’elles 


soient publiées on trouvera x et y d'après les tables bien connues donnant v en fonction de M. 


Représentation du lieu moyen. 


Pour l’époque 1925 Avril 8.1039 on trouve t-Tt--3.0779 jours, M= 0.032195, m = 0.032184, 
x—1.10468, y — 0.07120, et, d'après la formule (12), x= + 0.49773, y= — 0.88370, z =+ 0.44356. On a 
encore x + X=1.45148, y + Y= — 0.60338, z + Z= + 0.56516, d’où a= + 0.86894, b= P e 
c = + 0.33834. En comparant ces valeurs aux &, b © données, on obtient 

0-C = 0.00000 + 0.00001 — 0.00002. S 


Comme le calcul a été fait avec 5 décimales (excepté ses parties plus délicates), le lieu moyen 
est representé exactement, dans les limites des erreurs inévitables du calcul. 


Ephéméride. 

D’après les éléments ci-dessus on a léphéméride suivante (coordonnées vraies 1925.0 ; on obtien- 
drait facilement l’éphéméride rapportée à l’équinoxe du jour, en se servant de la Table de la correction 
des coordonnées rectangulaires du Supplemento Internationale Nr. 3 du Rocznik Astronomiczny, mais pour 
une éphéméride approchée les différences seraient insensibles). 


1925 Avril 16.0 @jgsso = 22738743 isz = +28°%50" r = 1.1203 p = 1.5708 


4 Ci 200 44 9 +33 51 1.1325 1.5282 
» 24.0 50 25 +39 6 1.1483 1.4927 
Le, 07 48 +44 36 1.1674 1.4654 
Mai 2.0 23 6 42 +50 14 1.1896 1.4476 


Les neuf constantes P, ... À, simplifient d’une manière essentielle la résolution du problème de 
la correction d’une orbite approchée, parce qu’elles permettent d'éviter l'emploi des constantes trigono- 
métriques de Schönfeld ou autres analogues. Ne voulant pas retarder l'apparition du présent article, 
nous remettons à une autre occasion l'exposition numérique de la méthode en question, reposant sur l'emploi de 
formules de nouveau genre. Le problème a été d’ailleurs déjà traité succinctement dans notre note „On 
a certain mathematical notion and its astronomical applications“ dans le Bull. de l Acad. Polon. d. Sciences, 
Série A. 1923. 


Cracovie, le 16 Avril 1925. Th. Banachiewicz. 


Contenu du Nr. 16: Th. Banachiewicz. La Comète Orkisz (1925 C), sa découverte, son orbite, formules de nouveau 
genre appliquées à la détermination d’une orbite. 


Drukowano na zlecenie Obserwatorjum Krakowskiego. 


1925. IV. 17. Dyrektor: Tad. Banachiewicz. 


